Gene expression analyses reveal differences in children's response to malaria according to their age.

In Bandiagara, Mali, children experience on average two clinical malaria episodes per year. However, even in the same transmission area, the number of uncomplicated symptomatic infections, and their parasitemia, can vary dramatically among children. We simultaneously characterize host and parasite gene expression profiles from 136 Malian children with symptomatic falciparum malaria and examine differences in the relative proportion of immune cells and parasite stages, as well as in gene expression, associated with infection and or patient characteristics. Parasitemia explains much of the variation in host and parasite gene expression, and infections with higher parasitemia display proportionally more neutrophils and fewer T cells, suggesting parasitemia-dependent neutrophil recruitment and/or T cell extravasation to secondary lymphoid organs. The child's age also strongly correlates with variations in gene expression: Plasmodium falciparum genes associated with age suggest that older children carry more male gametocytes, while variations in host gene expression indicate a stronger innate response in younger children and stronger adaptive response in older children. These analyses highlight the variability in host responses and parasite regulation during P. falciparum symptomatic infections and emphasize the importance of considering the children's age when studying and treating malaria infections.

Tick hemocytes have a pleiotropic role in microbial infection and arthropod fitness.

Uncovering the complexity of systems in non-model organisms is critical for understanding arthropod immunology. Prior efforts have mostly focused on Dipteran insects, which only account for a subset of existing arthropod species in nature. Here we use and develop advanced techniques to describe immune cells (hemocytes) from the clinically relevant tick Ixodes scapularis at a single-cell resolution. We observe molecular alterations in hemocytes upon feeding and infection with either the Lyme disease spirochete Borrelia burgdorferi or the rickettsial agent Anaplasma phagocytophilum. We reveal hemocyte clusters exhibiting defined signatures related to immunity, metabolism, and proliferation. Depletion of phagocytic hemocytes affects hemocytin and astakine levels, two I. scapularis hemocyte markers, impacting blood-feeding, molting behavior, and bacterial acquisition. Mechanistically, astakine alters hemocyte proliferation, whereas hemocytin affects the c-Jun N-terminal kinase (JNK) signaling pathway in I. scapularis. Altogether, we discover a role for tick hemocytes in immunophysiology and provide a valuable resource for comparative biology in arthropods.

Single-cell analyses of polyclonal Plasmodium vivax infections and their consequences on parasite transmission.

Most Plasmodium vivax infections contain genetically distinct parasites, but the consequences of this polyclonality on the development of asexual parasites, their sexual differentiation, and their transmission remain unknown. We describe infections of Saimiri monkeys with two strains of P. vivax and the analyses of 117,350 parasites characterized by single cell RNA sequencing and individually genotyped. In our model, consecutive inoculations fail to establish polyclonal infections. By contrast, simultaneous inoculations of two strains lead to sustained polyclonal infections, although without detectable differences in parasite regulation or sexual commitment. Analyses of sporozoites dissected from mosquitoes fed on coinfected monkeys show that all genotypes are successfully transmitted to mosquitoes. However, after sporozoite inoculation, not all genotypes contribute to the subsequent blood infections, highlighting an important bottleneck during pre-erythrocytic development. Overall, these studies provide new insights on the mechanisms regulating the establishment of polyclonal P. vivax infections and their consequences for disease transmission.

Gene expression analyses reveal differences in children's response to malaria according to their age.

In Bandiagara, Mali, children experience on average two clinical malaria episodes per season. However, even in the same transmission area, the number of uncomplicated symptomatic infections, and their parasitemia, vary dramatically among children. To examine the factors contributing to these variations, we simultaneously characterized the host and parasite gene expression profiles from 136 children with symptomatic falciparum malaria and analyzed the expression of 9,205 human and 2,484 Plasmodium genes. We used gene expression deconvolution to estimate the relative proportion of immune cells and parasite stages in each sample and to adjust the differential gene expression analyses. Parasitemia explained much of the variation in both host and parasite gene expression and revealed that infections with higher parasitemia had more neutrophils and fewer T cells, suggesting parasitemia-dependent neutrophil recruitment and/or T cell extravasation to secondary lymphoid organs. The child's age was also strongly correlated with gene expression variations. Plasmodium falciparum genes associated with age suggested that older children carried more male gametocytes, while host genes associated with age indicated a stronger innate response (through TLR and NLR signaling) in younger children and stronger adaptive immunity (through TCR and BCR signaling) in older children. These analyses highlight the variability in host responses and parasite regulation during P. falciparum symptomatic infections and emphasize the importance of considering the children's age when studying and treating malaria infections.

Gene expression analyses reveal differences in children's response to malaria according to their age.

In Bandiagara, Mali, children experience on average two clinical malaria episodes per season. However, even in the same transmission area, the number of uncomplicated symptomatic infections, and their parasitemia, vary dramatically among children. To examine the factors contributing to these variations, we simultaneously characterized the host and parasite gene expression profiles from 136 children with symptomatic falciparum malaria and analyzed the expression of 9,205 human and 2,484 Plasmodium genes. We used gene expression deconvolution to estimate the relative proportion of immune cells and parasite stages in each sample and to adjust the differential gene expression analyses. Parasitemia explained much of the variation in both host and parasite gene expression and revealed that infections with higher parasitemia had more neutrophils and fewer T cells, suggesting parasitemia-dependent neutrophil recruitment and/or T cell extravasation to secondary lymphoid organs. The child's age was also strongly correlated with gene expression variations. Plasmodium falciparum genes associated with age suggested that older children carried more male gametocytes, while host genes associated with age indicated a stronger innate response (through TLR and NLR signaling) in younger children and stronger adaptive immunity (through TCR and BCR signaling) in older children. These analyses highlight the variability in host responses and parasite regulation during P. falciparum symptomatic infections and emphasize the importance of considering the children's age when studying and treating malaria infections.

Tick extracellular vesicles impair epidermal homeostasis through immune-epithelial networks during hematophagy.

Hematophagous ectoparasites, such as ticks, rely on impaired wound healing for skin attachment and blood feeding. Wound healing has been extensively studied through the lens of inflammatory disorders and cancer, but limited attention has been given to arthropod-borne diseases. Here, we used orthogonal approaches combining single-cell RNA sequencing (scRNAseq), flow cytometry, murine genetics, and intravital microscopy to demonstrate how tick extracellular vesicles (EVs) disrupt networks involved in tissue repair. Impairment of EVs through silencing of the SNARE protein vamp33 negatively impacted ectoparasite feeding and survival in three medically relevant tick species, including Ixodes scapularis. Furthermore, I. scapularis EVs affected epidermal γδ T cell frequencies and co-receptor expression, which are essential for keratinocyte function. ScRNAseq analysis of the skin epidermis in wildtype animals exposed to vamp33-deficient ticks revealed a unique cluster of keratinocytes with an overrepresentation of pathways connected to wound healing. This biological circuit was further implicated in arthropod fitness when tick EVs inhibited epithelial proliferation through the disruption of phosphoinositide 3-kinase activity and keratinocyte growth factor levels. Collectively, we uncovered a tick-targeted impairment of tissue repair via the resident γδ T cell-keratinocyte axis, which contributes to ectoparasite feeding.

Tick hemocytes have pleiotropic roles in microbial infection and arthropod fitness.

Uncovering the complexity of systems in non-model organisms is critical for understanding arthropod immunology. Prior efforts have mostly focused on Dipteran insects, which only account for a subset of existing arthropod species in nature. Here, we describe immune cells or hemocytes from the clinically relevant tick Ixodes scapularis using bulk and single cell RNA sequencing combined with depletion via clodronate liposomes, RNA interference, Clustered Regularly Interspaced Short Palindromic Repeats activation (CRISPRa) and RNA-fluorescence in situ hybridization (FISH). We observe molecular alterations in hemocytes upon tick infestation of mammals and infection with either the Lyme disease spirochete Borrelia burgdorferi or the rickettsial agent Anaplasma phagocytophilum. We predict distinct hemocyte lineages and reveal clusters exhibiting defined signatures for immunity, metabolism, and proliferation during hematophagy. Furthermore, we perform a mechanistic characterization of two I. scapularis hemocyte markers: hemocytin and astakine. Depletion of phagocytic hemocytes affects hemocytin and astakine levels, which impacts blood feeding and molting behavior of ticks. Hemocytin specifically affects the c-Jun N-terminal kinase (JNK) signaling pathway, whereas astakine alters hemocyte proliferation in I. scapularis. Altogether, we uncover the heterogeneity and pleiotropic roles of hemocytes in ticks and provide a valuable resource for comparative biology in arthropods.

Malian children infected with Plasmodium ovale and Plasmodium falciparum display very similar gene expression profiles.

Plasmodium parasites caused 241 million cases of malaria and over 600,000 deaths in 2020. Both P. falciparum and P. ovale are endemic to Mali and cause clinical malaria, with P. falciparum infections typically being more severe. Here, we sequenced RNA from nine pediatric blood samples collected during infections with either P. falciparum or P. ovale, and characterized the host and parasite gene expression profiles. We found that human gene expression varies more between individuals than according to the parasite species causing the infection, while parasite gene expression profiles cluster by species. Additionally, we characterized DNA polymorphisms of the parasites directly from the RNA-seq reads and found comparable levels of genetic diversity in both species, despite dramatic differences in prevalence. Our results provide unique insights into host-pathogen interactions during malaria infections and their variations according to the infecting Plasmodium species, which will be critical to develop better elimination strategies against all human Plasmodium parasites.

Long read single cell RNA sequencing reveals the isoform diversity of Plasmodium vivax transcripts.

Plasmodium vivax infections often consist of heterogenous populations of parasites at different developmental stages and with distinct transcriptional profiles, which complicates gene expression analyses. The advent of single cell RNA sequencing (scRNA-seq) enabled disentangling this complexity and has provided robust and stage-specific characterization of Plasmodium gene expression. However, scRNA-seq information is typically derived from the end of each mRNA molecule (usually the 3'-end) and therefore fails to capture the diversity in transcript isoforms documented in bulk RNA-seq data. Here, we describe the sequencing of scRNA-seq libraries using Pacific Biosciences (PacBio) chemistry to characterize full-length Plasmodium vivax transcripts from single cell parasites. Our results show that many P. vivax genes are transcribed into multiple isoforms, primarily through variations in untranslated region (UTR) length or splicing, and that the expression of many isoforms is developmentally regulated. Our findings demonstrate that long read sequencing can be used to characterize mRNA molecules at the single cell level and provides an additional resource to better understand the regulation of gene expression throughout the Plasmodium life cycle.

Determination of the Stage Composition of Plasmodium Infections from Bulk Gene Expression Data.

Malaria symptoms are caused by the development of the parasites within the blood of an infected host. Bulk RNA sequencing (RNA-seq) of infected blood can reveal interactions between parasites and the host immune system during an infection, but because multiple developmental stages with distinct transcriptional profiles are concurrently present in infected blood, it is necessary to correct such analyses for differences in cell composition among samples. Gene expression deconvolution is a statistical approach that has been developed for inferring the cell composition of complex tissues characterized by bulk RNA-seq using gene expression profiles from reference cell types. Here, we describe the evaluation of a species-agnostic reference data set that can be used for efficient and accurate gene expression deconvolution of bulk RNA-seq data generated from any Plasmodium species and for correct gene expression analyses for biases caused by differences in stage composition among samples. IMPORTANCE Differences in cell type proportions among samples can introduce artifacts in gene expression analyses and mask genuine differences in gene regulation. Gene expression deconvolution allows estimation of the proportion of each cell type present in one sample directly from bulk RNA sequencing data, but this approach requires a reference data set with the signature profile of each cell type. Here, we evaluate the suitability of a rodent malaria parasite gene expression data set for estimating the proportions of each parasite developmental stage present in bulk RNA sequencing data generated from blood-stage infections with the human parasites Plasmodium falciparum and Plasmodium vivax. These analyses provide a species-agnostic approach for reliably estimating stage proportions in infected human blood and correcting subsequent gene expression analyses for these variations.

STRIDE: a command-line HMM-based identifier and sub-classifier of Plasmodium falciparum RIFIN and STEVOR variant surface antigen families.

BACKGROUND: RIFINs and STEVORs are variant surface antigens expressed by P. falciparum that play roles in severe malaria pathogenesis and immune evasion. These two highly diverse multigene families feature multiple paralogs, making their classification challenging using traditional bioinformatic methods. RESULTS: STRIDE (STevor and RIfin iDEntifier) is an HMM-based, command-line program that automates the identification and classification of RIFIN and STEVOR protein sequences in the malaria parasite Plasmodium falciparum. STRIDE is more sensitive in detecting RIFINs and STEVORs than available PFAM and TIGRFAM tools and reports RIFIN subtypes and the number of sequences with a FHEYDER amino acid motif, which has been associated with severe malaria pathogenesis. CONCLUSIONS: STRIDE will be beneficial to malaria research groups analyzing genome sequences and transcripts of clinical field isolates, providing insight into parasite biology and virulence.

Successful Profiling of Plasmodium falciparum var Gene Expression in Clinical Samples via a Custom Capture Array.

var genes encode Plasmodium falciparum erythrocyte membrane protein-1 (PfEMP1) antigens. These highly diverse antigens are displayed on the surface of infected erythrocytes and play a critical role in immune evasion and sequestration of infected erythrocytes. Studies of var expression using non-leukocyte-depleted blood are challenging because of the predominance of host genetic material and lack of conserved var segments. Our goal was to enrich for parasite RNA, allowing de novo assembly of var genes and detection of expressed novel variants. We used two overall approaches: (i) enriching for total mRNA in the sequencing library preparations and (ii) enriching for parasite RNA with a custom capture array based on Roche's SeqCap EZ enrichment system. The capture array was designed with probes based on the whole 3D7 reference genome and an additional >4,000 full-length var gene sequences from other P. falciparum strains. We tested each method on RNA samples from Malian children with severe or uncomplicated malaria infections. All reads mapping to the human genome were removed, the remaining reads were assembled de novo into transcripts, and from these, var-like transcripts were identified and annotated. The capture array produced the longest maximum length and largest numbers of var gene transcripts in each sample, particularly in samples with low parasitemia. Identifying the most-expressed var gene sequences in whole-blood clinical samples without the need for extensive processing or generating sample-specific reference genome data is critical for understanding the role of PfEMP1s in malaria pathogenesis. IMPORTANCE Malaria parasites display antigens on the surface of infected red blood cells in the human host that facilitate attachment to blood vessels, contributing to the severity of infection. These antigens are highly variable, allowing the parasite to evade the immune system. Identifying these expressed antigens is critical to understanding the development of severe malarial disease. However, clinical samples contain limited amounts of parasite genetic material, a challenge for sequencing efforts further compounded by the extreme diversity of the parasite surface antigens. We present a method that enriches for these antigen sequences in clinical samples using a custom capture array, requiring minimal processing in the field. While our results are focused on the malaria parasite Plasmodium falciparum, this approach has broad applicability to other highly diverse antigens from other parasites and pathogens such as those that cause giardiasis and leishmaniasis.

High-throughput detection of eukaryotic parasites and arboviruses in mosquitoes.

Vector-borne pathogens cause many human infectious diseases and are responsible for high mortality and morbidity throughout the world. They can also cause livestock epidemics with dramatic social and economic consequences. Due to its high costs, vector-borne disease surveillance is often limited to current threats, and the investigation of emerging pathogens typically occurs after the reports of clinical cases. Here, we use high-throughput sequencing to detect and identify a wide range of parasites and viruses carried by mosquitoes from Cambodia, Guinea, Mali and the USA. We apply this approach to individual Anopheles mosquitoes as well as pools of mosquitoes captured in traps; and compare the outcomes of this assay when applied to DNA or RNA. We identified known human and animal pathogens and mosquito parasites belonging to a wide range of taxa, as well as DNA sequences from previously uncharacterized organisms. Our results also revealed that analysis of the content of an entire trap could be an efficient approach to monitor and identify rare vector-borne pathogens in large surveillance studies. Overall, we describe a high-throughput and easy-to-customize assay to screen for a wide range of pathogens and efficiently complement current vector-borne disease surveillance approaches.

Variation in selective constraints along the Plasmodium life cycle.

Plasmodium parasites, the cause of malaria, have a complex life cycle, infecting alternatively vertebrate hosts and female Anopheles mosquitoes and undergoing intra- and extra-cellular development in several organs of these hosts. Most of the ~5000 protein-coding genes present in Plasmodium genomes are only expressed at specific life stages, and different genes might therefore be subject to different selective pressures depending on the biological activity of the parasite and its microenvironment at this point in development. Here, we estimate the selective constraints on the protein-coding sequences of all annotated genes of rodent and primate Plasmodium parasites and, using data from scRNA-seq experiments spanning many developmental stages, analyze their variation with regard to when these genes are expressed in the parasite life cycle. Our study reveals extensive variation in selective constraints throughout the parasites' development and highlights stages that are evolving more rapidly than others. These findings provide novel insights into the biology of these parasites and could provide important information to develop better treatment strategies or vaccines against these medically-important organisms.

Whole-genome analysis of Malawian Plasmodium falciparum isolates identifies possible targets of allele-specific immunity to clinical malaria.

Individuals acquire immunity to clinical malaria after repeated Plasmodium falciparum infections. Immunity to disease is thought to reflect the acquisition of a repertoire of responses to multiple alleles in diverse parasite antigens. In previous studies, we identified polymorphic sites within individual antigens that are associated with parasite immune evasion by examining antigen allele dynamics in individuals followed longitudinally. Here we expand this approach by analyzing genome-wide polymorphisms using whole genome sequence data from 140 parasite isolates representing malaria cases from a longitudinal study in Malawi and identify 25 genes that encode possible targets of naturally acquired immunity that should be validated immunologically and further characterized for their potential as vaccine candidates.

Transcriptional heterogeneity and tightly regulated changes in gene expression during Plasmodium berghei sporozoite development.

Despite the critical role of Plasmodium sporozoites in malaria transmission, we still know little about the mechanisms underlying their development in mosquitoes. Here, we use single-cell RNA sequencing to characterize the gene expression profiles of 16,038 Plasmodium berghei sporozoites isolated throughout their development from midgut oocysts to salivary glands, and from forced salivation experiments. Our results reveal a succession of tightly regulated changes in gene expression occurring during the maturation of sporozoites and highlight candidate genes that could play important roles in oocyst egress, sporozoite motility, and the mechanisms underlying the invasion of mosquito salivary glands and mammalian hepatocytes. In addition, the single-cell data reveal extensive transcriptional heterogeneity among parasites isolated from the same anatomical site, suggesting that Plasmodium development in mosquitoes is asynchronous and regulated by intrinsic as well as environmental factors. Finally, our analyses show a decrease in transcriptional activity preceding the translational repression observed in mature sporozoites and associated with their quiescent state in salivary glands, followed by a rapid reactivation of the transcriptional machinery immediately upon salivation.

Primaquine for Plasmodium vivax radical cure: What we do not know and why it matters.

Plasmodium vivax radical cure requires the administration of a blood schizonticide for killing blood-stage parasites and the addition of a drug able to kill hypnozoites, the dormant parasite stages residing in the liver of infected patients. All drugs used clinically for killing hypnozoites are 8-aminoquinolines and among them, primaquine has been at the forefront of P. vivax case management for decades. We discuss here the possible factors that could lead to the emergence and selection of P. vivax primaquine resistant parasites and emphasize on how a better understanding of the mechanisms underlying primaquine treatment and hypnozoite biology is needed to prevent this catastrophic scenario from happening.

Relative contributions of various endogenous and exogenous factors to the mosquito microbiota.

BACKGROUND: The commensal microbiota of mosquitoes impacts their development, immunity, and competency, and could provide a target for alternative entomological control approaches. However, despite the importance of the mosquito/microbiota interactions, little is known about the relative contribution of endogenous and exogenous factors in shaping the bacterial communities of mosquitoes. METHODS: We used a high-throughput sequencing-based assay to characterize the bacterial composition and diversity of 665 individual field-caught mosquitoes, as well as their species, genotype at an insecticide resistance locus, blood-meal composition, and the eukaryotic parasites and viruses they carry. We then used these data to rigorously estimate the individual effect of each parameter on the bacterial diversity as well as the relative contribution of each parameter to the microbial composition. RESULTS: Overall, multivariate analyses did not reveal any significant contribution of the mosquito species, insecticide resistance, or blood meal to the bacterial composition of the mosquitoes surveyed, and infection with parasites and viruses only contributed very marginally. The main driver of the bacterial diversity was the location at which each mosquito was collected, which explained roughly 20% of the variance observed. CONCLUSIONS: This analysis shows that when confounding factors are taken into account, the site at which the mosquitoes are collected is the main driver of the bacterial diversity of wild-caught mosquitoes, although further studies will be needed to determine which specific components of the local environment affect bacterial composition.

Host and Parasite Transcriptomic Changes upon Successive Plasmodium falciparum Infections in Early Childhood.

Children are highly susceptible to clinical malaria, and in regions where malaria is endemic, their immune systems must face successive encounters with Plasmodium falciparum parasites before they develop immunity, first against severe disease and later against uncomplicated malaria. Understanding cellular and molecular interactions between host and parasites during an infection could provide insights into the processes underlying this gradual acquisition of immunity, as well as to how parasites adapt to infect hosts that are successively more malaria experienced. Here, we describe methods to analyze the host and parasite gene expression profiles generated simultaneously from blood samples collected from five consecutive symptomatic P. falciparum infections in three Malian children. We show that the data generated enable statistical assessment of the proportions of (i) each white blood cell subset and (ii) the parasite developmental stages, as well as investigations of host-parasite gene coexpression. We also use the sequences generated to analyze allelic variations in transcribed regions and determine the complexity of each infection. While limited by the modest sample size, our analyses suggest that host gene expression profiles primarily clustered by individual, while the parasite gene expression profiles seemed to differentiate early from late infections. Overall, this study provides a solid framework to examine the mechanisms underlying acquisition of immunity to malaria infections using whole-blood transcriptome sequencing (RNA-seq).IMPORTANCE We show that dual RNA-seq from patient blood samples allows characterization of host/parasite interactions during malaria infections and can provide a solid framework to study the acquisition of antimalarial immunity, as well as the adaptations of P. falciparum to malaria-experienced hosts.

Single-cell transcription analysis of Plasmodium vivax blood-stage parasites identifies stage- and species-specific profiles of expression.

Plasmodium vivax and P. falciparum, the parasites responsible for most human malaria worldwide, exhibit striking biological differences, which have important clinical consequences. Unfortunately, P. vivax, unlike P. falciparum, cannot be cultivated continuously in vitro, which limits our understanding of its biology and, consequently, our ability to effectively control vivax malaria. Here, we describe single-cell gene expression profiles of 9,215 P. vivax parasites from bloodstream infections of Aotus and Saimiri monkeys. Our results show that transcription of most P. vivax genes occurs during short periods of the intraerythrocytic cycle and that this pattern of gene expression is conserved in other Plasmodium species. However, we also identify a strikingly high proportion of species-specific transcripts in late schizonts, possibly associated with the specificity of erythrocyte invasion. Our findings provide new and robust markers of blood-stage parasites, including some that are specific to the elusive P. vivax male gametocytes, and will be useful for analyzing gene expression data from laboratory and field samples.